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The  Validation  and  Application  of  a 
Bistatic  Two-Scale  of  Surface 
Roughness  Scattering  Model 


1.  INTRODUCTION 

Several  new  radar  systems  are  being  considered  which  operate  in  a  blstatic  mode.  This  offers  the 
potential  for  increased  target  radar  cross  sections  but  has  a  more  complicated  clutter  condition  than 
typical  monostatic  operation  would  involve.  This  report  addresses  three  topics  related  to  this 
question.  First,  a  bistatic  terrain  scattering  model  that  includes  two  distinct  sets  of  scattering 
phenomena1  is  evaluated  by  comparisons  with  some  bistatic  scattering  data.2  Secondly,  the 
contributions  from  each  of  the  levels  of  scattering  are  analyzed  as  a  function  of  bistatic  configuration. 
Finally,  the  validated  bistatic  scattering  results  are  compared  to  a  corresponding  set  of  monostatic 
values  representing  a  clutter  map  of  a  region  where  the  surface  is  loam. 


2.  MODEL  VALIDATION 

In  this  section  the  two-scale-of-roughness  bistatic  scattering  model  will  be  described.  Then  the 
data  base  used  for  the  comparisons  will  be  discussed.  Finally,  the  results  for  data  and  theory  will  be 


(Received  for  Publication  25  October  1988) 

1.  Ruck,  G.T.,  Barrick,  D.E.,  Stuart,  W.D.,  and  Krichbaum,  C.K.  (1970)  Radar  Cross  Section 

Handbook,  Vol.  2,  Plenum  Press,  New  York. 

2.  Cost,  S.T.  (1965)  Measurements  of  the  Bistatic  Echo  Area  of  Terrain  at  X-band,  Master’s  Thesis, 

Ohio  State  University. 
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examined  as  a  function  of  the  bistatic  scattering  angles.  Figure  1  shows  the  angles  involved.  9g  is  the 
complement  of  0j.  the  Incident  elevation  angle;  0S  is  the  scattered  elevation  angle;  and  4»s  is  the 
azimuthal  scattering  angle.  In  the  monostatic  case  4>s  =  180°. 

2.1  Theory 

The  model  is  quite  complex.  Details  are  given  here  to  make  the  comparisons  with  data  clearer  in 
terms  of  the  use  of  parameter  values  associated  with  the  particular  scattering  processes.  Two  distinct 
roughness  elements  are  assumed  to  contribute  to  the  scattering  cross  section  ct°  used  in  the  scattering 
model.  The  large  surface  roughness  is  described  by  a  model  based  upon  physical  optics  (PO) 
assumptions.  It  can  be  shown  that  one  sufficient  condition  for  the  validity  of  physical  optics  is  that 
TL  »  X,  where  TL  is  the  surface  correlation  length  associated  with  the  large  scale  roughness.  The  form 
for  a°  is  in  terms  of  the  Rayleigh  parameter  Z  =  y  (cos  0(  +  cos  9S).  Here.  oL  =  standard  deviation  in 
surface  height  for  large  scale  roughness.  When  the  conditions  for  the  validity  of  physical  optics  are 
met.  and  the  Rayleigh  parameter  Z  »1,  the  form  of  ct°  is  equivalent  to  a  high  frequency,  geometrical 
optics  (GO)  limit  solution.  The  small  scale  surface  roughness  is  described  by  the  perturbation  method 
(PM)  solution. 1  The  conditions  for  the  validity  of  the  perturbation  solution  are  that  the  rms  surface 
height  be  small  compared  to  a  wavelength  (2 n/X)  cts  <  1  and  that  the  surface  slopes  be  small  as/Ts  <  1 , 
where  as  =  standard  deviation  in  surface  height  for  the  small  scale  roughness  and  Ts  =  surface 
correlation  length  for  small  scale  roughness.  In  general,  the  surface  can  be  quite  complex;  its 
scattering  contribution  may  have  to  be  integrated  over  segments  with  different  geological  features. 
Here,  we  are  concerned  only  with  the  behavior  of  a°.  The  surface  areas  are  assumed  to  have  Gaussian 
height  statistics.  The  surface  correlation  function  is  assumed  to  be  Gaussian.  Surface  shadowing 
effects  are  also  included  in  the  formalism. 

For  physical  optics  models.  Ruck  et  al1  give  expressions  for  the  average  bistatic  rough  surface 
cross  section  cr°  under  the  following  assumptions:  (1)  the  radius  of  curvature  of  the  surface 
irregu laxities  is  larger  than  a  wavelength;  (2)  the  roughness  is  isotropic  in  both  surface  dimensions;  (3) 
die  correlation  length  is  smaller  than  either  the  x  or  y  dimension  of  the  sample  subregion;  and  (4) 
mtiltiDle  scattering  is  neglected.  Using  their  notation,  one  finds  that  the  expression  for  a°  becomes 

=  I  Ppq  I2js  (D 

where  flpq  represents  the  scattering  matrix  contributions,  S  is  the  local  shadowing  function,  and  the 
lerrn  J  is  related  to  the  surface  height  distributions  and  the  surface  slopes.  The  shadowing  function 
clearly  depends  on  the  roughness  of  the  surface,  and  Introducing  this  factor  into  the  analysis  can  have 
significant  effects  on  the  diffuse  scattered  power. 

In  this  report,  the  physical  optics  part  of  the  theoretical  cross  section  o°  is  not  evaluated  in  the 
high  frequency  geometrical  optics  limit.  Instead,  it  is  evaluated  by  making  the  assumption  of  small 
surface  slopes  at  J TL  <  1  so  that  J  is  expressed  as  a  single  integral  (Papa  et  al3). 


3  Papa,  R.J..  Lennon,  J.F.,  and  Taylor,  R.L.  (1984)  An  Analysis  of  Physical  Optics  Models  for  Rough 
Surface  Scattering.  RADC-TR-84- 195,  ADA  154960. 
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SURVEILLANCE  GEOMETRIES 


SATELLITE 


Figure  1.  Surveillance  Geometries  for  a  Bistatic  System  or  a  Monostatic  System 


3 


|  J0(vxyT>  lx2-  X  i  X  *  ! x  dt  (2) 

o 

where 

vxy  =  Vvx2  +  vy2 

vx  =  (2;t  /X) 

Vy  =  (271  /X) 

i  y2 

X  t  -  e  1  .  univariate  characteristic  function  of  the  surface  height  distribution  function 

X2  =  exp  |-I2(  1  -  e  XT  )]  bivariate  characteristic  function  (Gaussian  in  this  report) 

.Jn  is  the  zero  order  Bessel  function. 

i fere. 

qx  =  sin  0[  -  sin  9S  cos  <t>s  . 

4V  =  sin  0S  sin  0S  , 

=  -cos  0,  -  cos  0S  . 

It  is  in  the  scattering  matrix  term  that  the  dielectric  constant  representing  the  respective 
moisture  content  levels  is  Introduced.  In  that  term,  the  matrix  elements  for  linear  polarization  states 
are 

Pw  =  [ fi 2*i 3 R  ( 0 j " )  +  sin  0j  sin0s  sin2  <t> s  Ri(0|")l/[a  j a 4 ]  , 

Phv  =  sin  0S  |-sin  0(a3R  (0,")  +  sin  0S  a2RjJ0|”)l/la ia4)  .  (3) 

Pvh  =  sin  (sin  0sa2R  l0i")  _  sin  9t  a 3R_l (®i ” )  1  /( a  j a4)  , 

Pun  =  I- sin  0,  sin  0S  sin20s  R  ,  ( 0 j " )  -  a2a3Rj_(0(" )I/[a  j a4|  . 

Here,  R  ( 0 ( “ )  andRj_(0|")  are  Fresnel  reflection  coefTfcients 

R  (0,“)  =  [ercos0,-  -  yJtT[ir  -  sin2  0("  )/(  £rcos0j"  -  yjtT\ir  -  si  n2  Oj"  ]  (1) 

.10(1 

R  i  (0t " )  =  (  prcos  Oj"  -  *y/crhr  -  sin2  0)”  \/\  prcos  0("  +  'Jt^r  -  sin2  0,"  ]  . 
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Note  that  er  is  the  relative  complex  dielectric  constant  of  the  surface,  the  subscript  n  refers  to  the 
E-field  in  the  plane  of  incidence,  and  the  subscript  1  refers  to  the  E-field  normal  to  the  plane  of 
incidence.  The  remaining  angle-related  terms  are 

cos  0j"  =  (1  /V2" )  yj  1  -  sin  0j  sin  0S  cos  <{is  +  cos  0j  cos  0S 

a,  =  1  +  sin  0t  sin  0S  cos  <)>s  -  cos  0t  cos  0S  . 

a2  =  cos  ©!  sin  0S  +  sin  0j  cos  0S  cos  <}> s  . 


a3  =  sin  0j  cos  0S  +  cos  0t  sin  0S  cos  4>s  , 

and 

a4  =  cos  0j  +  cos  0S  . 

For  small  scale  roughness  the  cross  section,  o°ss  is  obtained  by  a  perturbation  method  solution 
to  Maxwell’s  equations.  The  fundamental  assumptions  for  this  case  are  small  roughness  (kos  <  1)  and 
small  surface  slopes  (I  3zs /3x  I .  I  dzs/dy  l<  1)  with  isotropic  roughness.  Here,  k  =  2 n/\,  and  Zs  is  the 
surface  height. 

For  this  model  we  have 

<r°ss  =  (4/n)k4  a|  cos2  0,  cos2  0S  I  apq  |2I  (5) 

where  Opq  is  proportional  to  the  scattering  matrix  element  and.  for  a  Gaussian  surface  with  Gaussian 
correlation 


The  a 

aim  = 


I  =  7tT2  exp  [-0.25k2Tf  (£  2  +  ^y2)]  . 
pq  terms  are  given  by 

l(pr  -  l)(pr  sin  0j  sin  0S  -cos  <j>s  yje.r\ir  -  sin2  0j  /erpr  -  sin2  0S  )  +  p2r(er  - 1)  cos  <)>s  ] 
[  pr  cos  0j  +  yjtT\ir  -  sin2  0j  ]  [  pr  cos  0S  +  -\/erpr  -  sin2  0S  ] 


(6) 


(7) 


aVII  -  sin<(»s 


er(pr  -  1)  Verpr  -  sin2  0t  -  pr(er  -  1)  Verpr  -  sin20s 
[pr  cos  0j  +  -\/erpr  -  sin2  0j  ]  (er  cos  09  +  yjt rpr  -  sin2  0S  1 


er(pr  -  1)  yj£r\ir  -  sin2  09  -  pr(er  -  1)  ^ernr  -  sin2  0, 

auv  =  si  n<j>9  , - - — .  , -  , 

[  er  cos  0(  +  Vert1r  -  sln20|  J  (  Hr  cos  +  VerHr  ~  sin2  09  J 


and 


avv 


[  (er  -  1) (e r  sin  0t  sin  09  -cos  <|>s  -\/er p.r  -  sin2  0t  -^erHr  -  sin2  09  )  +  e2r(pr  -  1)  cos  $s] 
[  er  cos  0j  +  Verllr  _  sln2  ®i  I  l  er  cos  0S  +  VerHr  -  sln2  9s  1 
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In  the  previous  equations.  er  is  the  relative  permittivity  of  the  surface  and  p,.  is  the  relative 
permeability.  They  may  be  either  real  or  complex  (for  a  lossy  surface).  Here  we  take  p,.  =  1. 

We  have  derived  individual  expressions  for  the  cross  section  of  the  surface  of  each  of  the  two 
levels  of  roughness  to  be  considered.  It  has  been  shown  that,  as  long  as  aL/TL  <  1  for  the  large  scale 
roughness,  the  total  scattering  cross  section  is  Just  the  sum  of  the  two  components.  If  that  is  not  the 
ease.  Brown4  has  shown  that  the  calculation  of  the  composite  cross  section  is  no  longer  that  simple. 
For  our  purposes  we  will  address  only  cases  that  do  satisfy  the  criterion,  aL/ Tl  <  1,  so  that  summation 
can  be  applied. 

In  this  report,  the  surface  heights  are  assumed  to  have  a  Gaussian  distribution.  Hence,  the 
shadowing  can  be  described  (in  the  high  frequency  limit  X  ->  0)  by  a  shadowing  function  S  derived  by 

Saucer.5 

(1)  For  <J)S  =  180° . 

S  =  1/(1  +  C0 )  when  9S  <  0t  ; 

S  =  1/(1  +  C2 )  when  0S  >  0t  ;  and  (8) 

(2)  For  <fs  *  180°. 

S  =  1/(1  +  C0+  C2)  . 

In  these  expressions, 

C0  =  [a,  tan  0t  /(T^)l  exp  {-[T/(2oL  tan  0,)12}  +  (1/2)  erfc  {t/(2ctl  tan  0,)} 

and 

C2  =  |aL  tan  0S  /(tVtc)]  exp  | - [T/(2ol  tan  8S)]2}  +  (1  /2)  erfc  |t/(2ol  tan  0S) J 

Brown4  has  shown  that  this  shadowing  function  S  multiplies  the  cross  section  for  both  large  scale 
roughness  and  the  cross  section  for  the  small  scale  roughness. 

2.2  Data  and  Model  Parameters 

Cost2  presents  blstatic  data  in  two  formats.  For  <{>s  =  0°,  values  of  o°  are  plotted  as  a  function  of 
elevation  scattering  angle  0S,  for  fixed  Incident  elevation  angles  0j.  In  the  second  series  of  data  sets  the 
values  of  cf  arc  plotted  as  a  function  of  azimuthal  scattering  angle  for  fixed  values  of  04  =  0S.  Both 
vertical- vertical  and  horizontal-horizontal  polarized  signals  were  considered  and  the  data  are  for  a 


4.  Brown,  G.S.  (1978)  Backscattering  from  a  Gaussian  distributed,  perfectly  conducting,  rough 
surface.  IEEE  Trans.  Antennas  Propag.  AP-26:472-482. 

5  Sancer,  M.I.  (1969)  Shadow  corrected  electromagnetic  scattering  from  a  randomly  rough  surface, 
IEEE  Trans.  Antennas  Propag.  AP-17:577-585. 
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wavelength  X  =  3  cm.  (X-band).  Only  limited  ground  truth  was  given  for  the  data.  For  the  azimuthal 
variations  we  will  show  comparisons  for  loam.  Cost  gives  the  standard  deviation  in  surface  heights 
for  his  loam  data  to  be  on  the  order  of  a  wavelength.  The  4>s  =  0°,  elevation  angle  cr°  comparisons  in 
this  report  are  for  sand.  Cost  gives  the  standard  deviation  of  the  surface  heights  for  sand  as  several 
wavelengths.  For  the  loam  comparisons,  three  sets  of  elevation  angles  were  considered:  9,  =  0S  =  60°. 

0j  =  0S  =  70°.  and  0j  =  0S  =  85°.  The  <J>S  =  0°  sand  comparisons  are  made  for  two  sets  of  incident  angles. 

9i  =  20°  and  0S  =  40°. 

The  model  parameters  have  to  be  related  to  the  scattering  surfaces  that  resulted  in  the 
experimental  data  sets.  The  parameters  must  have  physically  realistic  values  and  are  required  to 
satisfy  a  series  of  constraints  imposed  by  the  model.  First,  the  surface  must  be  very  large  compared  to 
the  large  scale  correlation  length.  Next,  the  large  scale  surface  slopes  must  be  small,  so  that  TL  >  aL. 
The  small  scale  surface  parameters  must  satisfy  the  requirements  that  Ts  «  TL  and  as  «  oL.  Also,  the 
small  scale  surface  parameters  must  satisfy  two  additional  constraints:  (1)  ^  as  <  1  and 
as/Ts  <  1.  Subject  to  all  these  conditions,  some  variation  in  the  surface  is  possible  and  there  was  no 
additional  ground  truth  supplied  by  Cost.  The  principal  model  parameter  of  interest  is  the  correlation 
length  Tl.  The  constraints  for  this  parameter  for  loam  are  such  that  0.06  m  <  TL  <  0.6  m  and  for  sand 
0. 15  m  <  TL  <  0.6  m  .  The  final  selection  of  parameters  was  made  by  considering  several  possibilities 
and  checking  agreement  with  a  limited  data  set.  The  same  parameters  were  then  used  to  predict 
behavior  for  other  cases  and  the  results  are  assessed  in  the  next  section.  As  an  example,  for  the 
azimuthal  cases,  small  incident  elevation  angle  cases  with  horizontally  polarized  signals  were  used  to 
establish  the  basic  parameter  set  consistent  with  the  parameter  constraints  of  the  model.  For  loam 
the  final  parameter  values  are  ctl2  =  0.03  m  (given);  TL  =  0.15  m  (mid-range);  as2  =  0.0045  m;  Ts  =  0.025 
m ;  and  e  =  30.0  +  J  2.0.  For  sand,  the  corresponding  baseline  parameters  are  slightly  different.  Here. 
ctl2  =  0.09  m  (given);  TL  =  0.45  m  (mid-range);  as2  =  0.0001  m;  Ts  =  0.0005  m;  and  e  =  3.0  +  J  0.2. 

2.3  Comparisons  of  Model  Results  with  Data 

The  two-scale-of-roughness  model  predictions  are  compared  to  the  X-band  data  taken  by  Cost. 
The  loam  data  allows  azimuthal  variations  to  be  examined  while  the  sand  results  show  scattering 
elevation  variations  for  <)>s=  0°. 

Figures  2  through  7  show  Cost’s  data  for  loam  as  points  for  0t  =  0S  =  60°,  70°.  and  85°  for  both  HH 
and  W  polarizations.  The  solid  lines  show  the  theoretical  model  results  with  the  two  scales  of 
roughness  parameters  chosen  by  judicious  variations  over  a  range  of  conditions.  The  dots  show  the 
experimental  data.  The  final  selection  of  parameters  has  been  cited  previously.  The  agreement 
between  the  experimental  data  and  the  theoretical  model  is,  in  general,  good,  considering  the  large 
variations  in  angles  of  incidence,  scattering  angles  and  polarizations.  It  should  be  noted  that  the 
general  trend  of  <r°  vs  <j>s  is  followed  by  both  the  experimental  data  and  the  theoretical  model.  For 
small  <t»s  (<J>S  «  0°),  o°  is  at  a  maximum;  as  4>s  increases,  o°  reaches  a  minimum  (20°  <  <f>s  <  100°);  and  then 
a°  increases  and  reaches  a  plateau  beyond  <|>s  =  100°.  This  is  true  for  0t  =  0S  =  60°  and  70°.  in  terms  of 
both  the  experimental  data  and  the  theoretical  model,  for  both  polarizations.  The  agreement  between 
the  experimental  data  and  the  theoretical  model  is  the  least  good  for  0,  =  0S  =  85°.  This  is  probably  due 
to  the  fact  that  here  multiple  scattering  plays  a  role,  and  the  two-scale  model  (with  shadowing)  doesn't 
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account  for  multiple  scattering.  Also,  the  nulls  in  o°  shown  In  the  theoretical  curves  near  <|>s  -  30°  In 
Figures  5  and  6  do  not  appear  in  the  experimental  data.  This  may  occur  because  the  nulls  are  narrow 
and  there  are  simply  not  enough  experimental  data  points  in  this  region. 

Figures  8  through  1 1  show  Cost’s  data  for  sand  as  points  for  0{  =  20°  and  40°  for  both  HH  and  W 
polarizations.  Again,  the  solid  lines  show  the  theoretical  model  results  for  the  two  scales  of  roughness 
parameters  chosen  in  accord  with  Cost's  data.  The  agreement  between  the  experimental  data  and  the 
theoretical  model  is  quite  good,  and,  in  general,  better  than  the  agreement  for  loam.  The  greatest 
discrepancy  between  the  experimental  data  and  the  theoretical  model  occurs  for  0S  ->  90°.  This 
discrepancy  again  may  be  attributed  to  the  effects  of  multiple  scattering,  which  mostly  occur  near 
small  grazing  angles.  Also,  there  may  be  large  errors  in  the  experimental  data  when  0S  -»  90°.  because 
it  is  difficult  to  define  precisely  the  radar  footprint  on  the  earth’s  surface. 

3.  ASSESSMENT  OF  SURFACE  ROUGHNESS  CONTRIBUTIONS 

In  the  previous  section  we  have  shown  how  the  two-scale  model  scattering  predictions  behave  as  a 
function  of  geometry.  In  this  section  we  will  examine  the  results  more  explicitly.  The  predictions  will 
be  examined  with  the  contribution  of  each  level  of  roughness  Isolated  from  the  other  so  that  the 
regions  where  one  or  the  other  scattering  mechanism  dominates  the  cross  section  can  be  seen  and 
analyzed.  Only  horizontal-horizontal  (HH)  polarization  results  will  be  considered,  since  the  behavior 
for  the  vertical- vertical  (W)  polarization  is  very  similar.  Figures  12  through  20  show  ct°  vs  <|>s  plots  for 
HH  polarization  for  the  large  scale  of  roughness  only  (loam  is  the  dielectric  surface).  Here,  the  surface 
slopes  are  small.  <tl/Tl  =  0.333  for  Figures  12.  13,  and  14;  and  cfl/Tl  =  0.577  for  Figures  15,  16,  and  17. 

The  figures  show,  in  general,  a  characteristic  trend.  The  normalized  cross  section  o°  starts  out  at 
a  relatively  high  value  at  <t>s  =  0°.  drops  down  to  some  minimum  value  as  4>s  increases,  and  then  rises  to 
a  plateau  as  <t>s  further  increases.  This  behavior  is  not  always  exhibited  when  0j  =  0S  =  85°,  where  o° 
often  exhibits  a  monotonic  decrease  as  <|>s  increases.  The  more  general  behavior  which  is  shown  at 
0,  =  0S  =  60°  and  70°  follows  what  Is  observed  experimentally.  Comparing  Figure  15  with  Figure  12, 
Figure  16  with  Figure  13,  and  Figure  17  with  Figure  14,  one  may  observe  another  general  trend  In 
behavior.  As  the  surface  slope  aL/TL  increases,  cr°  decreases  near  q>s  =  0°  and  o°  increases  for  <J>S  >  the 
null  angle. 

Figures  18.  19,  and  20  summarize  clearly  the  effects  of  oL  and  TL  on  the  graphs  of  or°  vs.  <|>s  for 
fixed  0j.  Comparing  Figure  19  with  Figure  18,  one  may  observe  a  general  trend:  as  TL  increases,  the 
right  hand  side  (backscatter)  decreases  and  the  left  hand  side  (forward  scatter)  increases  (for  ctl 
constant,  both  HH  and  W  polarizations).  Comparing  Figure  20  with  Figure  19,  we  may  note  that  as  oL 
decreases,  the  right  hand  side  (backscatter)  decreases  and  the  left  hand  side  (forward  scatter)  increases. 
Thus,  for  the  large  scale  roughness,  both  TL  and  aL  control  the  shape  of  the  o°  vs.  <(>s  curves. 

In  Figures  21  through  29,  only  the  small  scale  roughness,  (os,Ts)  contributes  to  the  normalized 
cross  section  o°.  Comparing  Figure  24  with  Figure  2 1 ,  Figure  25  with  Figure  22,  and  Figure  26  with 
Figure  23,  one  may  observe  a  trend  in  behavior  that  is  analogous  to  the  large-scale-of-roughness  case. 
As  Ts  decreases  (and  os/Ts,  the  slope,  increases),  a0  decreases  near  <t>9  =  0°  and  o°  increases  for  $s  >  30°. 
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Figure  10.  RCS  vs.  0S  for  Sand 
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Figure  18.  RCS  vs.  <(>s  for  Loam 
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Figures  27.  28,  and  29  summarize  the  effects  of  the  small  scale  roughness  on  the  o°  vs.  $s  curves. 
Comparing  Figure  28  with  Figure  27.  one  may  note  a  general  trend;  as  Ts  Increases,  the  right  hand  side 
(backscatter)  decreases  and  the  left  hand  side  (forward  scatter)  Increases  and  the  entire  angular  extent 
where  a°  contributes  becomes  narrower  (for  crs  constant,  both  HH  and  W  polarizations).  Thus,  Ts 
controls  the  shape  of  the  <r°  vs.  <)>s  curve.  Comparing  Figure  29  with  Figure  27,  one  may  note  that  as  os 
decreases,  the  entire  a°  vs.  <(>s  curve  moves  down  (decreases)  and  does  not  change  shape  (for  Ts  constant, 
both  HH  and  W  polarizations). 

4.  BISTATIC  AND  MONOSTATIC  CLUTTER  MAPS 

One  question  of  interest  Is  how  the  clutter  cross  section  map  of  a  region  would  change  if  a  bistatic 
map  were  considered  Instead  of  a  monostatic  one.  We  will  use  our  validated  bistatic  o°  model  to 
demonstrate  the  differences  for  a  particular  scattering  configuration.  The  region  will  consist  of  a 
loam  surface  extending  some  200  nmi  down  range  and  400  nmi  cross  range.  A  monostatic  system  is 
assumed  to  be  at  a  grazing  angle  of  about  6°  with  the  front  center  of  the  region  and  about  2°  with  the 
rear.  (The  grazing  angle  is  the  complement  of  9t).  These  conditions  apply  to  the  scattered  elevation 
angle  values  of  the  bistatic  system.  The  bistatic  Incident  elevation  angle  Is  taken  with  respect  to  the 
center  of  the  region,  0j  =  70°.  For  both  cases  the  cr°  maps  are  constructed  by  dividing  the  region  Into  5 
nmi  by  5  nmi  boxes  and  assigning  a  o°  value  to  the  center  of  each  box.  Since  the  region  is  so  extensive, 
the  further  refinement  of  a  spherical  geometry  was  Introduced.  The  four-thirds  earth  radius  concept 
then  takes  into  account  tropospheric  refraction  effects. 

For  this  section  the  frequency  Is  considered  to  be  L-band  (f  =  1  GHz).  The  terrain  is  considered  to 
be  very  rough  but  the  constraint  ctl/Tl  <  1  Is  maintained.  Here  aL  =  1.0  m:  TL  =  1. 1  m;  a8  =  0.0045  m, 

Ts  =  1  m;  and  e  =  30.0  +  J  2.0. 

In  Figure  30,  the  bistatic  clutter  map  for  loam  is  shown  for  horizontal-horizontal  polarization, 
with  the  wavelength  equal  to  30  cm.  The  first  observation  to  be  noted  is  that  the  clutter  cross  section 
a°  is,  in  general,  greatest  In  the  direct  forward,  near  specular,  direction.  This  is  due  to  the  scattering  In 
the  forward,  near  specular  direction  (small  slopes,  ctl/Tl  <  1).  The  decrease  in  clutter  cross  section  o° 
with  increasing  range  Is  due  primarily  to  shadowing.  The  surface  irregularities  cause  less  of  the  rough 
surface  to  be  illuminated  by  the  Incident  rays  as  the  elevation  scattering  angle  0S  decreases. 

In  Figure  31,  the  monostatic  clutter  map  for  loam  is  shown  for  horizontal-horizontal 
polarization.  It  may  be  noted  immediately  that,  for  a  given  clutter  cell,  the  monostatic  cross  section  a° 
is.  In  general,  at  least  50  dB  less  than  the  corresponding  bistatic  cross  section.  Since  the  monostatic 
cross  section  o°  has  a  fixed  azimuthal  scattering  angle  (<]>s  =  180°),  the  monostatic  clutter  map  Is 
circularly  symmetric  about  the  radar  position.  The  behavior  of  the  monostatic  cross  section  o°  with 
range  is  similar  to  the  bistatic  cross  section;  o°  decreases  as  the  range  increases  primarily  because  of 
shadowing. 

Figure  32  is  not  a  clutter  map,  rather,  it  is  a  clutter  ratio  map,  that  is,  o°HH  (monostatic) /a°HH 
(bistatic).  It  may  be  noted  Immediately  from  Figure  32  that  for  any  cell,  the  monostatic  cross  section 
is  at  least  50  dB  less  than  the  corresponding  bistatic  cross  section  <r°. 

Figure  33  is  a  bistatic  clutter  map  for  loam  for  vertical-vertical  polarization.  Comparing 
Figure  33  with  Figure  30  shows  that  in  the  forward  scatter  direction,  the  bistatic  cross  section  c^w  Is 
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Figure  31.  Monostatic  Clutter  Map  for  Loam  (HH  Polarization) 
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Figure  32.  Monostatic /Bistatic  Ratio  Map  for  Loam  (HH  Polarization) 
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about  20  dB  less  than  o0,^ .  For  larger  azimuthal  scattering  angles,  corresponding  to  a  cross  range 
>  50  nmi,  the  blstatic  cross  section  o° w  is,  in  general,  greater  than  o°HH.  The  behavior  of  o° w  with 
range  Is  similar  to  o°HH ;  it  decreases  as  the  range  Increases  because  of  shadowing. 

Figure  34  is  a  monostatic  clutter  map  for  loam  for  vertical-vertical  polarization.  Comparing 
Figure  34  with  Figure  31,  we  may  note  that  the  monostatic  cross  section  o°w  for  a  given  clutter  cell  is 
very  nearly  Identical  tc  cr°HH. 

Figure  35  is  a  map  of  the  ratio  <j°w  (monostatic) /a°w  (blstatic).  It  may  be  noted  immediately  that 
the  monostatic  cross  section  a°w  is  at  least  45  dB  less  than  the  blstatic  o°w  for  any  given  clutter  cell. 
Comparing  Figure  35  with  Figure  32,  It  may  be  noted  that  the  two  map  ratios  <r°(monostatlc)/ 
o°(bistatic)  for  the  two  polarizations  are  quite  distinct.  This  is  due  primarily  to  the  differences  in  the 
two  bistatic  cross  sections. 


5.  CONCLUSIONS 

Reasonable  agreement  was  obtained  between  the  blstatic  scattering  model  with  two  scales  of 
roughness  and  the  X-band  data.2  Both  elevation  plane  comparisons  for  sand  and  azimuthal  results 
for  loam  were  considered.  Considering  the  sparseness  of  the  data  and  the  limited  nature  of  the  ground 
truth,  it  appears  that  sufficient  accuracy  and  trend  preservation  can  be  obtained  using  the  theoretical 
model.  It  should  be  noted  that,  where  the  agreement  is  poorest  (small  incident  grazing  angle),  both  the 
data  and  the  model  are  suspect.  The  uncertainty  in  the  data  is  caused  by  the  footprint  determination 
requirement  while  the  physical  optics  model  does  not  allow  multiple  scattering  which  is  more  likely 
to  occur  at  those  conditions. 

The  parametric  studies  in  which  each  level  of  roughness  was  treated  separately  allowed  us  to 
evaluate  the  model  performance.  The  physical  optics  terms  were  dominant  In  the  forward  scattering 
directions  and  in  those  cases  the  results  were  controlled  by  the  surface  slope.  As  the  scatterirq^ 
extended  beyond  these  regions,  the  physical  optics  results  were  not  as  significant  and  the  behavior  was 
no  longer  slope  dependent.  Different  patterns  occurred  depending  on  whether  the  standard  deviation 
in  heights  or  the  correlation  length  was  varied.  The  perturbation  theory  (small  scale)  results  do  not 
show  any  direct  slope  related  pattern.  For  all  cases  there  is  a  distinct  difference  in  behavior  to  the 
curves  when  either  as  or  Ts  is  varied.  as  affects  the  magnitude  and  Ts  affects  the  shape  of  the  small 
scale  pattern. 

When  the  validated  bistatic  scattering  model  was  used  to  compare  blstatic  and  monostatic  clutter 
cross  sections  for  a  loam  covered  surface  subdivided  Into  5  nmi  boxes  the  results  showed  that  the 
bistatic  a°  values  for  the  configuration  used  in  the  analysis  always  exceed  the  corresponding 
monostatic  values  by  50  dB  for  horizontal  polarization  and  by  45  dB  for  vertical  polarization. 
Shadowing,  as  expected,  played  a  more  significant  role  for  the  furthest  down  range  cells.  For  this  type 
of  geometry,  any  enhanced  target  detection  using  bistatlcs  has  to  come  from  greatly  increased  target 
cross  sections  at  blstatic  angles. 
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Figure  34.  Monostatic  Clutter  Map  for  Loam  (W  Polarization) 
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